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I I Abstract. We present a study of small carbon chains and rings in Photon Dominated Regions (PDRs) performed at 

^ . millimetre wavelengths. Our sample consists of the Horsehead nebula (B33), the pOph L1688 cloud interface, and 

' the cometary-shaped cloud IC63. Using the IRAM 30-m telescope, the SEST and the Effelsberg 100-m telescope 

, at Effelsberg., we mapped the emission of C2H, C-C3H2 and C4H, and searched for heavy hydrocarbons such as 

CO ' C-C3H, /-C3H, Z-C3H2, /-C4H2 and CeH. The large scale maps show that small hydrocarbons are present until the 

' edge of all PDRs, which is surprising as they are expected to be easily destroyed by UV radiation. Their spatial 

, distribution reasonably agrees with the aromatic emission mapped in mid-IR wavelength bands. C2H and C-C3H2 

JnI ■ correlate remarkably well, a trend already reported in the diffuse ISM (Lucas and Liszt 2000). Their abundances 

I relative to H2 are relatively high and comparable to the ones derived in dark clouds such as L134N or TMC-1, 

. known as efficient carbon factories. The heavier species are however only detected in the Horsehead nebula at 

CL(' a position coincident with the aromatic emission peak around 7/im. In particular, we report the first detection 

JL I of CeH in a PDR. We have run steady-state PDR models using several gas-phase chemical networks (UMIST95 

■ and the New Standard Model) and conclude that both networks fail in reproducing the high abundances of some 
' of these hydrocarbons by an order of magnitude. The high abundance of hydrocarbons in the PDR may suggest 
, that the photo-erosion of UV-irradiated large carbonaceous compounds could efficiently feed the ISM with small 

IL" ' carbon clusters or molecules. This new production mechanism of carbon chains and rings could overcome their 

destruction by the UV radiation field. Dedicated theoretical and laboratory measurements are required in order 
, to understand and implement these additional chemical routes. 

H : 
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1. Introduction 

Carbon is the fourth most abundant element in the in- 
terstellar medium (ISM), and also the most versatile for 
building molecules. Carbon chemistry can therefore be 
considered as the core of interstellar chemistry. Of the 
nearly 130 molecules now observed in various sources, 
about 75% have at least one carbon atom, while one 
fourth are hydrocarbons. Moreover, the heaviest and most 
complex molecules are organic molecules with carbon. 
This statistic does not take into account the Polycyclic 
Aromatic Hydrocarbons (PAHs), nor the Diffuse Infrared 
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Band carriers (DIB) which are most likely large organic 
molecules (Herbig 1995). Carbon, neutral or ionized, is 
also one of the main reactants in interstellar chemistry 
networks, due the large number of organic molecules, but 
also to its versatility allowing it to participate to numer- 
ous chemical reactions at any temperature, from the very 
cold dense cores, to warm and hot gas. Therefore under- 
standing the carbon chemistry is of major importance for 
astrochemistry, and for star formation. 

A large fraction of the chemical reactions in ISM net- 
works involve the numerous hydrocarbons present in the 
ISM. Such molecules were first reported in circumstel- 
lar shells, where the chemistry is particularly favorable 
to their formation (C4H - Guelin et al. 1978, C5H - 
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Cernicharo et al. 1986), but also in molecular dark clouds, 
where they are efficiently shielded from the interstellar ra- 
diation (C3H2 - Thaddeus 1981 et al., C2H- Gottheb et al. 
1983). Since then, hydrocarbons of increasingly greater 
number of carbons have been reported in several objects 
(e.g. CgH - Guehn et al. 1987, CgH - Dickens et al. 2001). 
In the diffuse gas, molecules like CH, CH+ and CN are 
known since the 40's, and carbon clusters (C2, C3) are 
now almost routinely detected in the visible towards bright 
stars (Maier et al. 2001, Roueff et al. 2002, Oka et al. 
2003) as well as in the far-IR (e.g. Cernicharo et al. 2000). 
Using radio telescopes, Lucas & Liszt (2000) have also 
shown that C2H and C-C3H2 are ubiquitous in diffuse gas, 
confirming previous work by Cox et al. (1989). 

A natural question therefore arises: if carbon chains 
are present in the diffuse ISM, what does happen in the 
Photon-Dominated Regions (hereafter PDRs) ? As in dif- 
fuse clouds, the chemical processes are dominated by the 
radiation but the gas is denser. The radiation field can also 
be constrained (intensity and direction) and, with the new 
data obtained by ISO, an accurate picture of the mid-IR 
emission due to the Aromatic Infrared Band (AIB) car- 
riers at these cloud interfaces has emerged. Knowledge 
about the distribution of carbon chain and rings in PDRs 
is however yet scarce and limited to the works of Fuentc et 
al. (1993, 2003), who reported the observation of C2H at 
some positions in NGC7023, the Orion Bar and NGC7027, 
and Ungerechts et al. (1997), who presented maps of the 
Orion Bar in the millimetre transitions of C2H and cyclic- 
C3H2. In this paper we extend the study of hydrocarbons 
in PDRs through an extensive inventory of carbon chains 
and rings of up to six carbon atoms. Our aims are to 
get a better view of their distribution on large scale with 
respect to other standard tracers, and to derive their con- 
tribution to the total carbon budget in comparison with 
values measured in other regions. 

The paper is organised as follows. We present in Sect.|21 
the three sources studied here, while the observations are 
described in Sect.|31 We then analyse the spatial distribu- 
tion of the mapped hydrocarbons (Sect. I4.1|l and derive 
the molecular column densities of the observed species 
(Sect. I5.1|l . These results are compiled in the form of a 
carbon budget and compared to similar observations re- 
ported in dark clouds and in the diffuse ISM (Sects. 15.21 
and 16.2(1 . We finally compare the inferred abundances to 
numerical PDR models fScct. lT^ and present our conclu- 
sions in Sect.[7| 

2. Presentation of the sources 

2.1. The Horsehead nebula 

In the visible, the Horsehead nebula (B33, Barnard 1919) 
appears as a dark patch seen in silhouette against the 
ionized Ha emission emanating from the IC434 Hii re- 
gion. The large-scale maps reported by e.g. Maddalena 
et al. (1986, C0(2-l)) or Lada et al. (1991, CS) show 
its connection to the L1630 cloud, which belongs to the 
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Fig. 1. Composite colour image of the Horsehead nebula 
PDR obtained by the VLT in the B, V and R bands. The 
arrow indicates the direction of the illuminating star cOri. 
Image courtesy of ESO. 
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Fig. 2. Visible image from the Digital Sky Survey in the 
area of IC63. The cometary shape originates from the 
strong radiation field emitted by the close-by HD5394 
stars 20' South- West. 
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Fig. 4. Integrated emission maps (T^, Kkm/s) of the molecular emission in the Horsehead nebula PDR for several 
species. Contours are: c-C3H2(2i,2-lo,i) 0.1 to 0.7 by 0.15, C2H(l-0) 0.3 to 2.1 by 0.3, C4H(9-8) 0.025 to 0.2 by 0.03, 
C^®0(2-1) 1 to 6 by 1. Also displayed is the continuum emission between 5 and 8.5 /im (aromatic emission) as mapped 
by ISO (from Abergel et al. 2002, contours are 10 to 25 by 5 MJy/sr), and at 1.2mm (contours are 3 to 53 by 10 
mJy/11" beam) mapped by the MAMBO 117-channel bolometer. The white squares indicated the position labelled 
as CO-peak (top) and IR-peak (bottom). 



Orion B molecular complex (see e.g. Abergel et al. 2002 
for an overview of the region). In the scenario proposed 
by Rcipurth & Bouchet (1984), the Horsehead nebula is 
interpreted as an early Bok globule emerging from its 
parent cloud via the eroding incident radiation field emit- 
ted by the close-by 09.5V star uOri. The Horsehead neb- 
ula thus corresponds to a condensation resisting the il- 
lumination and it presents on its western border a PDR 
whose detailed shape is displayed on Fig. ^ Assuming a 
distance of 400 pc (Anthony-Twarog 1982), and consider- 
ing a projected distance to trOri of 0.5°, the incident radi- 
ation field is of order Go = 100 (Zhou et al. 1993, Abergel 
et al. 2003). This can be seen as relatively weak compared 
to other well-known PDRs of the solar neighborhood (e.g. 
the Orion Bar, NGC7023). This PDR was mapped by ISO 
(6" resolution) in filters around 7 and 15 //m (Abergel et al. 
2002, Abergel et al. 2003). It reveals a very thin filament 
(10", or '--^0.02pc at 400 pc) whose narrow size suggests 
that the PDR is seen perfectly edge-on (see also Fig.QJ. 



Little attention was paid to this object until recently. 
In the CO and isotope maps of Kramer et al. (1996), the 
nebula is heavily diluted in the 2' beam and only ap- 
pears as a small intrusion of LI 630 into the Hii region. 
These large-scale data already indicated mean densities of 
order 3 x lO^cm"'^ (Zhou et al. 1993) and column den- 
sities compatibles with Ay 3 (Kramer et al. 1996). 
Observations at higher spatial resolution (10-20") were 
only reported very recently. Pound et al. (2003) presented 
BIMA CO(1~0) and discussed in details the possible ori- 
gin and evolution of the dark cloud. They reveal a com- 
plex velocity structure but their data are likely affected by 
non- negligible optical depth effects. Abergel et al. (2003) 
report IRAM 30-m maps obtained in the (J=l-0) and 
( J=2-l) transitions of CO, ^^CO and C^^O. They infer 
densities compatible with the values reported at larger 
scale, and show clear selective photo-dissociation of the 
considered species at the PDR edge. 
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Fig. 3. Close-up of the western border of L1688 mapped 
by ISOCAM around 7 urn (Abergel et al. 1996). The 
illuminating star HD 147889 direction is indicated by 
an arrow. The black squares indicate the positions of 
the cut where molecular data have been collected. The 
tilted rectangle delimitates the area mapped in ^'^CO and 
Ci^O. Offsets are given with respect to a2ooo = 16''25'"58", 
^2000 - -24°21'00". 



2.2. IC63 

IC63 is a reflexion nebula associated to the BO.SIVpe vari- 
able star 7 Cas (HD5394, Go = 1100) located 230±70pc 
from us (Vakili et al. 1984). This strong radiation field 
has created a PDR at the western border of the cometary- 
shaped molecular cloud reported by Jansen et al. (1994). 
IC63 is associated to another reflexion nebula, namely 
IC59, located at a projected distance of 20' (or 1.3pc), 
as is illustrated on Fig. [3 According to Blouin et al. 
(1997), IC59 would be lying in the background of the IC63- 
HD5394 pair. Jansen et al. (1994), on the other hand, 
claim that the illuminating star would be closer to the 
observer, so that IC63 is partially seen face-on. 

The molecular component appears elongated and 
seems well aligned with the direction of 7 Cas. Jansen 
et al. (1994) report that the CO and CS emissions are 
confined into a cometary cloud of size '^l'x2', and as- 
sume an even smaller area of 30" x 20" for all other species 
than CO. This has particular consequences in terms of 
beam dilution. These authors derive Tj^in — 50 ± lOK and 
= 5 ± 2 X 10** cm^'^, with a possible density gradi- 
ent along the cloud major axis. The inferred H2 column 
densities translate into Ay = 6.3 ± 2.5 mag. Finally, an 
analysis of the carbon budget in this source indicates that 
the total gas phase carbon abundance (5.4xl0~^) is only 
13% of the solar one, which suggests that the bulk of the 



carbon is in solid phase and might be located on PAHs or 
larger grains. 

IC63 has also been mapped with ISOCAM in the 
CVF mode allowing imaging spectroscopy between 5 and 
16.5 /xm (see Fosse et al. 2000). Fig. shows the emis- 
sion observed around the 6.2 /im AIB feature, which cor- 
responds to a stretching mode of the C-C bond. The aro- 
matic emission follows the PDR border and exhibit several 
peaks just behind the PDR and inside the molecular tail. 

2.3. The p Ophiuchi L1688 cloud 

The Ophiuchi molecular complex is a close-by star for- 
mation region (distance 135±15pc, Habart et al. 2003). 
While its so-called "North" area is considered inactive, 
the West counterpart contains two main massive cores, 
L1688 and L1689, heated on large scales by the Sco-OB2 
association. We concentrate here on the L1688 complex, 
which hosts numerous pre-stellar cores (e.g. Motte et al. 
1998), as well as hundreds of YSOs (Greene & Young 
1992, Bontemps et al. 2001), very likely responsible for 
the strong luminosity in the infrared. 

ISOCAM observations revealed bright interlaced fila- 
ments of width ~ 0.03 pc coinciding with the western bor- 
der of the cloud (Fig. see also Abergel et al. 1996). It 
corresponds to a PDR (hereafter called L1688-W) illumi- 
nated by the B2V star IID147889. According to Abergel et 
al. (1999), this star is located in the centre of a spherical 
cavity exhibiting irregular edges. The filamentary struc- 
ture is interpreted as the edge-on regions of this cavity. 
Assuming that the star and PDRs are in a common plane 
perpendicular to the line of sight, Habart et al. (2003) in- 
fer an incident radiation field of Go — 400. Based on the 
AIB and the Il2(l-0) S(l) line emissions, these authors 
predict a density plateau of uh = 4 x lO'* cm~^ inside the 
cloud, followed by a decrease of the form nH(?') oc r^-^ to- 
wards the radiating star, where r is the radial dimension 
across the PDR. This drop occurs at Aa ~ —100" on the 
cut indicated in Fig. |21 

On large scales, the mid-IR emission shape is well re- 
produced by the CO and isotopes maps obtained by Lada 
& Wilking (1980) and Wilkingfe Lada (1983). Their data 
revealed a complex velocity structure and, in particular, 
significant self-absorption in the ^^CO species. On these 
scales, ^^CO data are compatible with Ay in the range 
50-100 mag. 

3. Observations 

The data presented in this paper have been gathered be- 
tween September 1999 and May 2002 at the IRAM 30-m, 
the CSO, the SEST and the Effelsberg lOO-m^ telescopes. 
They are concentrated on a series of molecules including 
CO and its main isotopomers, as well as a collection of 
small carbon chains and rings. The details on the con- 

^ The 100-m telescope is operated by the MPIfR (Max- 
Planck-Institut fiir Radioastronomie) 
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Molecule 


Transition 


Frequency 


HPBW 






(GHz) 


(arcsec) 


^''CO 


J=2-l 


230.538000 


11 


"CO 


J=2-l 


220.398686 


22/44 


ci«o 


J=l-0 


109.782160 


22/44 




J=2-l 


219.560357 


11 


CS 


J=2-l 


97.980950 


24 


C2H 


N=l-0, J=3/2-l/2 








F=2-l 


87.316925 


28/56 




F=l-0 


87.328624 


28/56 




N=l-0, J=l/2-l/2 








F=l-1 


87.402004 


28/56 




F=0-1 


87.407165 


28/56 




N=3-2, J=7/2-5/2 








F=4-3 


262.004260 


28('') 




F=3-2 


262.006480 


28(") 




N=3-2, J=5/2-3/2 








F=3-2 


262.064990 


28(") 




F=2-l 


262.067460 


28(") 


C-C3H 


N=2i2-lii, J=5/2-3/2 








F=3-2 


91.494349 


27/54 




F=2-l 


91.497608 


27/54 




N=2i2-lii, J=3/2-l/2 








F=l-0 


91.692752 


27/54 




F=2-l 


91.699471 


27/54 


^-CsH 


J=9/2-7/2 








F=5-4 (f) 


97.995166 


25/50 




F=4-3 (f) 


97.995913 


25/50 




F=5-4 (e) 


98.011611 


25/50 




F=4-3 (e) 


98.012524 


25/50 


C-C3H2 


li,o~lo,i 


18.343145 


54(6) 




2l,2~lo,l 


85.338898 


28/56 


^-C3H2 


1q \ — Oo 


20.792590 


48"'' 




33,0~22,1 


216.27875 


11 


C4H 


N=9-8, J=19'/2-17/2 








F=9-8 


85.634006 


28/56 




F=10-9 


85.634017 


28/56 




N=9-8, J=17/2-15/2 








F=8-7 


85.672581 


28/56 




F=8-8 


85.672583 


28/56 


/-C4H2 


lln 11 — l(Jn in 


98.245016 


24 


CH3C2H 


J=5-4 








K=2 


85.450730 


28 




K=l 


85.455622 


28 




K=0 


85.457272 


28 


CeH 


^113/2 J=59/2-57/2 e 


81.777893 


28 




^Hg/a J=59/2-57/2 f 


81.801247 


28 


HN"C 


J=l-0 








F=0-1 


87.090735 


27 




F=2-l 


87.090859 


27 




F=l-1 


87.090942 


27 


HCsN 


J=9-8 


81.881468 


28 


Table 1. 


Line parameters for the all the species and data 



reported in this paper. When two spatial resolutions are 
indicated, we refer to the HBPW's corresponding to the 
IRAM 30-m and the SEST telescopes respectively. (""^ : ob- 
served at the CSO. Observed at the EfFelsberg lOO-m 
telescope. See Sect.Olfor details. 



sidered lines are compiled in Tab. ^ We give here some 
additional information for each source. In the following 
the prefixes c- and I- will respectively refer to the cyclic 
and linear isomers of the considered species. 

3.1. Horsehead nebula 

The Horsehead nebula PDR was mapped at the 30- 
m telescope in the (J=1~0) and (J=2-l) transitions of 
C^^O, as well as in c-C3H2(2i,2-lo,i)> C2H(l-0) and 
C4H(9-8) (Fig. 0}. System temperatures at 3mm were 
in the range 100-120 K. Some weaker lines were probed 
at dedicated positions (Figs. |H1 and |5Jl in CS, c-C3H(2i_2- 
/-C3H(9/2~7/2), Z-C3H2(33,o-22,i), /-C4H2(llo,ii' - 
10o,io), HNi3c(l-0), CH3CCH(5k-4k) and C6H(59/2- 
57/2). Except from the on-the-fly maps of C^^O, all spec- 
tra were obtained in the Frequency Switching mode using 
an autocorrelator spectrometer providing 80 kHz resolu- 
tion elements. Additional data were obtained in C2H(3- 
2) at the CSO (AOS spectrometer, T^ys ^ BOOK), and in 
Z-C3H2(lo,i-Oo,o) and c-C3H2(li,o-lo,i) at the Effelsberg 
100-m telescope {Tgys ^ 90 K). Moreover, a 1.2 mm contin- 
uum emission map was obtained at the 30-m using the new 
117-channel MAMBO bolometer (11" resolution). Using 
a fast-mapping mode (Teyssier & Sievers 1999), we could 
cover areas of 8.5' x 7.5' in about 30 minutes, with differ- 
ent chopping throws and at different hour angles, reducing 
efficiently some mapping artifacts inherent to the EKH in- 
version technique applied here (Emerson et al. 1979). Only 
part of this map is displayed here (Fig. 0J and further re- 
sults from these data will be discussed in a forthcoming 
paper. 

The data were first calibrated to the scale using 
the so-called chopper wheel method (Penzias & Burrus 
1973). The final adopted scale however depends on the 
source size. In the particular case of the Horsehead nebula, 
we applied the correction factor introduced in Abergel et 
al. (2003) taking into account the error beam of the 30- 
m telescope. At 3 mm, we assumed hydrocarbon emission 
areas similar to the C^^O one. 

3.2. IC63 

IC63 was mapped at the 30-m in C0(2-l) and C2H(l-0) 
using the on-the-fly mode (Fig. ISJ . Orthogonal coverages 
were optimally combined using the PLAIT algorithm by 
Emerson & Graeve (1988). In addition to 0^^0(2-1), the 
same lines as in the Horsehead were probed along a cut 
crossing the PDR (see Fig. EJ. All other observing pa- 
rameters were similar to the Horsehead ones. Dedicated 
positions were also observed at the 100-m in c-C3H2(li,o~ 
lo,i)- 

After calibrating to the scale, the lack of large- 
scale information in the probed molecules motivated us to 
adopt the widely-used main beam brightness temperature 
Tmb assuming beam efficiencies at the IRAM 30-m of 1% 
(= Soff/Fcff in the IRAM nomenclature) = 0.79 and 0.49 
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at 3 mm and 1.3 mm respectively. An additional correction 
was considered to account for the source dilution in the 
beam (see the description by Jansen et al. 1994) applying 
deconvolved sizes as inferred in Sect. 14.31 

3.3. L1688-W 

The L1688-W cloud was observed at the SEST in July 
2000. We mapped the CO (2-1) and 01^0(1-0) line emis- 
sions in an area indicated on Fig.|31 The area we have cov- 
ered is just at the edge of the large-scale map of Wilking & 
Lada (1983). This figure also shows the position of the cut 
performed in Frequency Switching for some of the small 
hydrocarbon species listed above. Fig. |H1 displays the inte- 
grated intensity maps of the two CO isotopomers. The cor- 
responding observing parameters can be found in Tab. ^ 
As for 1C63, the data were finally scaled to main beam 
brightness temperature, assuming SEST beam efficiencies 
of 77b = 0.78 and 0.56 at 3 mm and 1.3 mm respectively. 

4. Results 

4.1. General trends 

Despite the strong incident radiation fields (Go ~ 100- 
1000) compared to the diffuse ISM (Go ~ 1), the small 
carbon chains and rings are detected in all objects. At 
first glance, this is surprising as small hydrocarbons are 
expected to be rapidly destroyed by such UV radiations. 
The spatial intensity variations of all hydrocarbons mea- 
sured in this study correlate well. This result was al- 
ready reported by Lucas & Liszt (2000) in the diffuse ISM 
for C-C3H2 and C2II. Fig. [7| illustrates this trend in the 
Horsehead nebula. Some additional information is given 
below for each of the sources. 

4.2. Horsehead nebula 

The maps obtained in the Horsehead nebula are displayed 
in Fig. 01 They show that the hydrocarbon large-scale dis- 
tribution correlates well with the PDR structure as traced 
by ISO or ^^CO (see maps by Abcrgcl et al. 2003). In par- 
ticular the sharp edge, although smoothed by the beams, 
is observed in all tracers. 

On smaller scales, different behaviours have to be men- 
tioned. It is indeed interesting to note the different po- 
sitions of the respective hydrocarbon and CO (and iso- 
topomers) peaks. Fig. [3 shows that the correlation be- 
tween C^^O and C2II is worse than with C-C3H2, though 
the former lines are the strongest and thus offer better 
signal-to-noise ratio. Moreover, the hydrocarbon peak is 
located very close to the peak of aromatic emission around 
7 /ini (hereafter called IR peak) , suggesting a chemical link 
between AIB carriers and the small hydrocarbons. More 
conclusive comparisons with the ISO mid-IR data are nev- 
ertheless limited by the 30-m spatial resolution. This has 
important consequences on the physico-chemical analysis 
along cuts crossing the PDR as, for instance, any abun- 
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Fig. 5. Maps of IC63. Top: Integrated emission of the 
ISOCAM CVF spectra between 5.98 and 6.54 ^m after 
subtraction of a 1-order baseline. Contours are 70 to 
310 by 40 MJy/sr. The dashed box indicates the area 
covered by the molecular millimetre data displayed be- 
low. Bottom: Integrated emission maps (T^, Kkm/s) of 
i2CO(2-l) (contours 5 to 40 by 5) and C2H(l-0) (con- 
tours 0.10 to 0.35 by 0.05). The white squares indicate 
the position of the IR peaks. 



dances probed at the IR peak will not coincide with the 
peak emission of the hydrocarbons. This issue will be ad- 
dressed in a companion paper presenting interferometric 
observations of the same molecules (resolution ~ 5", Pety 
et al. 2003). 

Other species were also probed at positions of inter- 
est. They include heavier and rarer hydrocarbons, as well 
as some cyanopolyynes or density tracers such as CS (see 
Tab.^ and were chosen in order to complete our inventory 
of the carbon budget. In the Horsehead, these deep inte- 
grations were unfortunately performed prior to the large- 
scale mapping, so that the position observed at the IR 
peak does actually not correspond to the maximum hy- 
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Fig. 6. Channel maps of the C^*O(l-0) emission (Kkm/s, 
TJ) in the L1688-W PDR, as well as integrated intensity 
maps of Ci**O(l-0) and i3C0(2-l) ((a) and (b)) in the 
velocity interval 3.4-4.2 km/s, corresponding to the emis- 
sion associated to the PDR as seen by ISO. The velocities 
of each channel map is indicated in the upper left cor- 
ner. Offsets are given with respect to a2ooo = 16^25™58'', 
<52ooo = -24°21'00". 



drocarbon emission, as is seen on Fig. 0] The CO peak 
position was inferred from the ^^CO intensity map ob- 
served by Abergel et al. (2003). Although not displayed in 
Fig- El HN^^C was detected at both positions, and HC3N 
was observed at the hydrocarbon peak. I-C3H2, ^-C4H2 
and CH3C2H were not detected at either of the positions, 
while C-C3H and ^-CsHwere detected only at the IR peak. 
All other molecules were detected to a better than So- 
level. We also emphasize the detection of CgH at the IR 
peak which is to our knowledge the first detection of this 
molecule reported in a PDR. It could however not be ob- 
served at the CO peak. 

4.3. IC63 

The IC63 maps are shown in Fig. [SJ as well as a cut il- 
lustrating the radial profiles of more species across the 
PDR border (Fig. ^J). The correlation observed between 
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Fig. 7. Scatter diagrams in the Horsehead nebula. The 
scatter around the zero emission points is indicative of 
the noise in the data. 



the various tracers in the Horsehead is not so striking in 
the IC63 maps as this source is barely resolved by the 30- 
m observations and the large beam (28" compared to the 
6" offered by ISO) significantly smoothes out the details 
at the PDR border. Smoothing the ISO data to the 30-m 
resolution however shows a very good agreement between 
the probed emissions as the two main infrared peaks get 
smoothed into one single maximum coincident with the 
molecular emission feature. On the other hand, the CO 
and C2II peaks coincide fairly well, and are located just 
North to the south-most aromatic peak at 6.2 /im (indi- 
cated on Fig. . From the observed integrated emission, 
a beam-deconvolved size of ^ 40"x30" (FWHM) can be 
inferred for ^^CO and C2H (the North-East tail is not 
accounted here). 

In IC63, the deep integrations were done along a cut 
crossing the south-most AIB emission peak (indicated in 
Fig. ISJ , as well as some extra positions inside the cloud 
tail. C4H and C-C3II2 are in most cases detected, the lat- 
ter including transitions at 18 and 85 GHz (see Tab. 
C2H(3~2) is only marginally detected at some positions. 
The non-detections are similar to the ones reported in 
the Horsehead nebula, but include as well both isomers of 
C3H, HNi^C and CgH. 

4.4. L1688-W 

In L1688-W, the molecular cloud exhibits a complex ve- 
locity field, as was already reported by Wilking & Lada 
(1983). It is illustrated on Fig. |Hlfor the area covered in 
this study. We associate the PDR border as seen by ISO 
(Abergel et al. 1999) to the molecular emission in the 
velocities range 3.4 to 4.2 km/s (see integrated intensity 
maps in Fig. EI). Along the cut displayed in Fig. |3 the 
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Fig. 10. Normalized line intensities across the PDR 
fronts of the Horsehead nebula and IC63. Horsehead: iso- 
declination cut performed through the IR-peak (see Fig. 2)) 
at ^2000 = — 2°28'00". All intensities are normalized to 
the peak (respectively 0.24, 0.9, 2.3 and 8.7 K from top 
to bottom, and 24MJy/sr for the ISO LW2 profile). The 
illuminating source lies towards the negative offsets. IC63: 
cut performed along the direction indicated by open cir- 
cles on Fig. El All intensities are normalized to the peak 
(respectively 0.03, 0.27, 0.5 and 0.2 K from top to bottom, 
and 310MJy/sr for the ISOCAM CVF profile). The illu- 
minating source lies towards the negative offsets. All error 
bars are ±1 a. Upper limits indicate non-detections at a 
2(7 level. 



brightest hydrocarbon emission is detected at positions 
towards the illuminating star where 0^*0 has already 
started to decay significantly (Fig. 111(1 . The emission at 
negative offsets can nevertheless be due to the particular 
geometry of L1688-W(Abergel et al. 1999), implying that 
the PDR is partly seen face-on at probed positions inside 
the cavity. 

In L1688-W, a somewhat smaller number of species 
was investigated in deeper integrations. Apart from the 
data presented in Fig. observations of C-C3H and l- 
C3H were unfruitful down a 5mK (Icr) noise level. 



5. Analysis of the physico-chemical conditions 

We have used the available collection of species to derive 
the physical and chemical conditions at work in the three 
sources, and in particular assess a tentative carbon budget 
in these PDRs (see Sect. l6.1|l . 

5.1. Column densities 

As mentioned in the previous section, one should bear in 
mind that beam dilution effects are expected to affect the 
measured intensity. This is especially true in our case since 
we look at positions located close to sharp edges. To first- 
order, this mainly means that some of the column densities 
derived below may be under-estimated. The assumptions 
and methods used to compute the quantities of interest 
are described in details in Appendix ^ We discuss in the 
following some specific hypothesis on a source to source 
basis. The results are gathered in Tab. El 

5.1.1. Horshead nebula 

From the ^^CO maps of Abergel et al. (2003), we derived 
Tkin 40 ± 2 K at the PDR location. Using this tempera- 
ture, the hydrogen number densities deduced from LVG 
modelling are of the order of some 10^ cm^"^. 

Following the method described in Appendix IA.2I we 
inferred C2II excitation temperatures between 5 and 10 K, 
while Ttot never exceeded 1.5. We cross-checked these re- 
sults using the additional N=3-2 transitions observed at 
some positions. The N=3-2/N=l-0 line intensity ratio are 
consistent with Tex in the range 4.5-8 K, in good agree- 
ment with the HFS outputs. The results show similar val- 
ues as the one reported by Fuente et al. (1993, 2003) 
in NGC7023 and the Orion Bar. The constant density 
assumed to conduct the C-C3H2 and C4II analysis (see 
Appendix I A. 3(1 was taken to be 2xl0'*cm~'^. 

5.1.2. IC63 

As discussed previously, the relatively small size of the 
emitting area compared to the mapping beams requires 
a correction factor to be applied to the line intensities. 
We used the approach described by Jansen et al. (1994) 
to convert T^b into brightness temperatures within a 
source filling 1000 arcsec^ for ^^CO and the hydrocar- 
bon emission, and 400 arcsec^ for the ^^C0(2-l) emission. 
Under these assumptions, the ^^CO observations indicate 
Tkin ^ 48±2K, in excellent agreement with Jansen et al.'s 
finding of 50 ± 10 K. Since only one C^^O transition was 
available in this source, we made the additional assump- 
tion of a volume density — 5±2x 10* cm~'^, based on 
the results of Jansen et al. (1994). The same density was 
applied to the analysis of the C-C3II2 and C4H column 
densities. The inferred C^^O column density at (0",0") 
offset is however 3 times lower than the value reported by 
these authors. This is likely due to the 10" offset between 
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Fig. 8. Sample spectra of heavier hydrocarbons and other species obtained at the so-called CO peak (indicated by 
dashed ^^CO contours, 30% and 80 to 95% of the peak, from Abergel et al. 2003) and IR peak (indicated by full line 
ISO contours 15 to 30MJy/sr, from Abergel et al. 2002). All temperature are in K, scale. The same scales are used 
for both positions to ease the comparison between the various molecules. The missing channels in the CgH spectrum 
were removed due to spurious features in the backend response. For this line, we also display the Gaussian fit to the 
line. Offsets are given with respect to a2ooo=05'^40'"59.0^ Jaooo = -02°25'31.4". 



their position and ours, sufficient to explain such rapid 
variations in regards of the source size (see Tab. 13). 

The C2H analysis indicates that 80% of the observed 
points exhibit Ttot < 1 and are consistent with Tex = 5.8- 
8 K. The inferred abundances are somewhat higher than 
the results inferred by Jansen et al. (1994) at a position 
10-20" away from ours. 

5.1.3. L1688-W 

The CO kinetic temperatures were derived using an ad- 
ditional C^^0(2-l) map obtained at the CSO (Gerin, 
private communication). We calculated Tkin ~ 21 ± 2K 
in the region of interest. As already reported by Lada 
& Wilking (1980), the CO lines however exhibit signif- 
icant self-absorption. Despite the dilution effects men- 
tioned above, the temperatures in this source are thus 
lower limits, implying upper limits for the molecular col- 
umn densities. 

While the C2H HFS analysis reveal excitation tem- 
perature similar to the ones found in IC63, we used the 
volume density profile derived by Habart et al. (2003) to 
constrain the C-C3H2 and C4H column densities. The sim- 



ulations confirm the small optical depth at the probed 
positions. 

5.2. Fractional abundances 
5.2.1. Systematics 

The molecular abundances with respect to H2 are com- 
piled in Tab. 12 Appendix IA.6l describes the methods used 
to derive the H2 column densities in each of the sources. 
The hydrocarbon abundances appear reasonably constant 
across the different objects and positions, although to a 
lesser extent in L1688-W where abundances do not always 
peak at the intensity maxima. We will discuss this point 
later. This constancy confirms the general trend already 
reported in the diffuse gas (Lucas & Liszt 2000). Although 
error bars remain high (essentially due to uncertainties on 
A^(H2)), the observed values result in average abundances 
of [C2H] = 1.9 ± 0.4 X 10-^ [C-C3H2] = 1.4 ± 0.2 X 10"^ 
and [C4H] = 1.1 ± 0.2 X 10"^. 

Abundance ratios between molecules show a somewhat 
similar trend and seem to be hierarchically distributed. On 
average, the following ratios are found: [C^^0]/[C2H] ^ 25 
(although with large scatter), [c-C3H2]/[C4H] ~ 1 and 
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Fig. 9. Sample spectra of weaker hydrocarbons and other species obtained at some dedicated positions in IC63 dis- 
played over contours of the 6.2 /.tm ISO integrated emission map. It illustrates in particular the non-detection of the 
cyclic and linear forms of C3H around 92 and 98 GHz respectively. The probed positions are marked with filled and 
open circles. The open circles indicate the cuts displayed in Fig. ^1 All temperature are in K, scale. Offsets are 
given with respect to Q^2ooo=00'^59"^00.7^ (52ooo = 60°53'19". 



[C2H]/[c-C3H2] ~ 15. The latter result is of the same order 
as the value reported in the diffuse gas (ratio of 28±8) by 
Lucas & Liszt (2000). 

5.2.2. Horsehead nebula 

In this source, it is interesting to note that the position 
where the hydrocarbon emission is the largest coincides 
with a peak of the H2 column density as traced by the con- 
tinuum emission at 1.2mm (Fig.Q. This trend does how- 
ever not hold for species such as C^^O (although within 
large error bars), CS and HN^'^C, whose abundance max- 
ima do not coincide with this II2 column density peak. 
Note also that the fraction of the beam filled by the 
cloud is smaller at the hydrocarbon peak than at the CO 
peak, thus lowering the intensities measured at the bor- 
der. This is confirmed by interferometric observations of 
the PDR border which show an enhancement of the hy- 
drocarbon abundances at the aromatic peak (Pety et al. 
2003). 

5.2.3. IC63 

In IC63, the small relative variations of abundances ob- 
served between the different offsets is very likely related 



to the proximity of the probed positions with respect to 
the beam sizes. The C^^O abundances are significantly 
lower than the value reported by Jansen et al. (1994). 
This discrepancy arises from the large II2 column densi- 
ties assumed at positions probed outwards the cloud. 

5.2.4. L1688-W 

Firm conclusions are more difficult to draw in L1688-W as 
the available sample remains restricted to few points. If 
our picture of the variations of the H2 column density 
is correct, the highest hydrocarbon abundances are not 
found at their corresponding intensity peaks but closer 
to the illuminating star (Fig. 111(1 . Again, this can be due 
to the particular geometry of the PDR partly seen face- 
on. Still, the C^^O abundances experience a drop much 
steeper than the hydrocarbons in the direction of the illu- 
minating star. 

6. Carbon chemistry 

6.1. Carbon budget 

The compilation of molecular abundances allows us to de- 
rive a tentative carbon budget representative of PDR con- 
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Fig. 11. Normalized line intensities over the cut displayed on Fig. O across the L1688-W PDR. Error bars are ±lcr. 
The star lies towards the negative offsets. The 0-offset refers to the most North-East point of the cut (Fig. O. All 
intensities are normalized to the peak (respectively 0.08, 0.2, 0.43 and 1.6K from top to bottom, and 117MJy/sr for 
the ISO LW2 profile). An upper limit is indicated for the rightmost point of C4H. Sample spectra at two offsets are 
also shown. The dashed line indicates the centroid velocity of the optically thin lines, illustrating the self-absorption 
in ^^CO. All temperature are in K, scale. 



ditions. As a matter of comparison. Tab. |31 also includes 
results obtained in the diffuse ISM (Lucas & Liszt 2000), 
as well as in a typical carbon-rich dark cloud (see also 
Sect. l6.2|l . For this, we have co-added typical abundances 
for molecules with a given number of carbon atoms. Since 
all molecules are not observed in the same sources, this 
carbon budget has only a statistical meaning, and num- 
bers cannot be very accurate. The contribution of car- 
bon clusters containing up to 2-3 atoms is very similar for 
all three media, with a clear abundance decrease as the 
molecular size increases, both in PDRs and dark clouds. 
Heavier chains and rings (4 carbon atoms and more) are 
however less abundant in PDRs and diffuse clouds than in 
dark clouds by a factor of order 5-10. The fraction of car- 
bon tied to more complex molecules may be only partly 
representative of the real situation. A lower number of 
molecules of that size is detected and other species still 
remain to be identified. 

6.2. Comparison to dark clouds 

We have compared the carbon chain and ring abundances 
to results obtained in TMC-1 (e.g. Pratap et al. 1997) and 



L134N (L183, e.g. Dickens et al. 2000), two dark clouds 
known as efficient carbon factories. The carbon budget 
corresponding to the dark cloud conditions is also indi- 
cated in Tab. |31 As for the PDRs, it remains difficult to 
derive accurate estimates of the total H2 column densities 
in these objects. These estimates are indeed usually af- 
fected by strong molecular depletion (e.g. Cernicharo et al. 
1987, Alves et al. 1999) and possible evolution of the grain 
properties in the cold and dense medium characteristic of 
dark clouds (e.g. Stepnik et al. 2003, Kramer et al. 2003). 
The abundances reported in Tab. |31 have been calculated 
assuming an averaged extinction of order 20 mag, which 
may lie in between the actual values applying to TMC- 
1 and L134N. In addition we have decided to compare 
the fractional abundances with respect to C2H, the most 
abundant hydrocarbon in our sample. Assuming again an 
extinction of 20 mag, the C2H fractional abundances to 
H2 at the TMC-1 cyanopolyyne peak (CP) and at the 
L134N hydrocarbon peak would however be 4 x 10~^ and 
8 X 10^^ respectively, somewhat lower than the values we 
have derived in PDRs. This is also true for species such 
as C-C3H2 or C4H. 
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Table 2. Column densities and abundances relative to H2 (boldface) inferred at positions of interest in the Horsehead, 
IC63 and L1688-W. Units for column densities are in cm^^ and indicated in superscript for each molecule. 

Units for abundances are indicated in the first column of each corresponding row. The offsets refer to the positions 
indicated Figs. El and El 



Tab. 0] compiles the abundances relative to C2H for 
each of the sources at the positions of interest. The error 
bars have not been considered as we look here at general 
trends. Contrary to the abundances relative to H2, the 
fractional abundances relative to C2H are generally larger 
in the dark clouds than in our sample PDRs (we will here 
not address the differences between TMC-1 and L134N, 
see Fosse 2003). For all molecules but Z-C3H, there is at 
least a factor of 5 difference between the molecular frac- 
tions in the PDRs and in TMC-1. This is particularly true 
for abundant species such as C-C3H2 and C4H, whereas 
they seem to have similar contributions to the carbon bud- 
get in both media (Tab. |3J). This behaviour is very likely 
due to the different chemistries at work in the two media, 
and will be further discussed in the next section. We also 
note the enormous drop of [HC3N]/[C2H] between TMC-1 
CP and the Horsehead nebula (factor of 1000). Although 
we probably missed the "cyanopolyyne" peak position of 
the Horsehead, the difference is still meaningful. 

6.3. Comparison to models 

In this section, we compare our observed data with model 
predictions. We used a plane parallel PDR model which 
solves self-consistently the thermal and chemical balance 
of a semi-infinite homogeneous slab (Le Bourlot ct al. 
1993, Le Petit et al. 2002). We describe the results for 
physical conditions mimicking those of the Horsehead neb- 



ula and the p Oph interface, namely a gaz density of njj = 
2 X IQ'^cm"'^ and a radiation field multiplied by Go — 100 
relative to the ISRF determined by Habing (1968). Since 
the chemical networks have not been fully tested for large 
molecules, we have chosen to use two different networks, 
the New Standard Model developed by E. Herbst and col- 
laborators (Terzieva & Herbst 1998) and the UMIST95 
chemical rate file (Millar et al. 1997). Both chemical 
networks include about 400 hundred species with 
up to 8 carbon atoms, and 4000 chemical reactions. 
The various possible isomers are not differenciated 
in the chemical networks. 

The photodissociation of H2 and CO and the 
photoionization of C are treated specifically in the 
PDR code via the computation of the actual ultra- 
violet radiative transfer in the presence of the com- 
plex chemistry provided by the two chemical net- 
works. Then, the transition of atomic to molecular 
hydrogen is identical in both models. However, the 
balance between neutral carbon and carbon con- 
taining species in the intermediate range of Ay is 
dependent both on the photo-processes (identical 
in the two networks) and the individual chemical 
reactions which involve different assumptions. 

Except for the formation of molecular hydrogen, we 
have not considered any reaction on grain surfaces. All 
molecules are therefore formed and destroyed in the gas 
phase. For the external layers of the cloud where the radi- 
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1.1 X 10" 


10 




4.2 X 10"i° 


7 (HC7N,..) 








5.5 X 10"^° 


8 (CgH,..) 








1.1 X 10"" 



Table 3. Carbon budget derived in PDRs (this work), diffuse ISM (from Bell et al. 1983, Federman et al. 1994, 
Shuping et al. 1999, Lucas & Liszt 2000, Roueff et al. 2002, Cernicharo et al. 2002 and Rachford et al. 2002) and dark 
clouds (from Pratap et al. 1997, Bell et al. 1998, Dickens et al. 2000, Turner et al. 2000, Fosse et al. 2001 and Fosse 
2003). Note that, for a given number of carbon atoms, the molecules found in the different regions may differ (e.g. 
no C„ carbon clusters in PDRs and dark clouds). For the difTuse ISM, we used the C/H abundance ratio of 
1.4x10"'* reported by Cardelli et al. (1996) and assumed an H/H2 abundance ratio of 0.6 as observed 
in the C Oph line of sight (e.g. Rachford et al. 2002). 



Molecule 


Horsehead 


IC63 


L1688-W 


TMC-1 


L134N 




IR peak area 


(0",-0") 


(Aq=-21") 


(CP) 




C2H (10^3 cm"2) 


16 


5.3 


3.0 


7.2 


15 


C-C3H 


0.024 


<0.023 




0.140 


0.050 


/-C3H 


0.013 


<0.110 




0.012 


0.040 f''^ 


C-C3H2 


0.058 


0.085 


0.100 


0.800 


0.330 


;-C3H2 


<0.002 






0.030 


0.004 (''^ 


CH3C2H 


<0.069 






1.4 


0.070 '•''^ 


C4H 


0.069 


0.062 


0.500 


2.8 


0.735 '''^ 


CeH 


0.004 






0.115 


0.008 '"^^ 


HN^^^C 


0.002 


<0.002 




0.056 


0.054 


HC3N 


0.003 






2.2 (-f) 


0.127 



Table 4. Fractional abundances relative to C2H. The assumed column densities of C2H are recalled in the first row. 
CP stands for Cyanopolyyne Peak. References are: ('^^ Pratap et al. 1997, ^''^ Fosse et al. 2001, Guehn et al. 1982, 
('^^ Fosse 2003, ('^^ Dickens et al. 2000, Takano et al. 1998. Since no CgH and HC3N measurements were 
available at the Horsehead IR peak, the corresponding abundances relative to C2H are given for the 
hydrocarbon peak. 



ation field is the strongest, the main destruction process 
is photo-dissociation, while chemical reactions are more 
important in well shielded regions. PAHs are taken into 
account for their role in the gas heating, but do not par- 
ticipate in the chemistry. Using two networks helps us 
in setting error bars on model predictions and therefore 
in identifying agreement or discrepancy between our ob- 
served data and the model predictions more securely. 

The PDR model assumes steady state equilibrium for 
both the thermal processes and the chemistry. Given the 
large density and radiation field, this assumption should 
be justified since the chemical time scales are dominated 
by the photo-dissociation reactions which are very rapid 
in the conditions we modelled. Furthermore, there is no 
need to invoke time dependent chemistry for the synthesis 
of molecules since ionised and atomic carbon are present 
with a large abundance in the external layers (up to Ay 



= 4 - 5 mag for C), and can participate in a large variety 
of chemical reactions for building molecules. 

We show in Fig. El model predictions for the two 
chemical networks separately. Though the transition from 
atomic to molecular hydrogen is very similar for both 
chemical networks, the transition from ionised carbon to 
neutral carbon, then to carbon monoxide, shows small dif- 
ferences, with the formation of carbon monoxide occurring 
at slightly smaller Ay in the NSM compared to UMIST95. 

We now discuss the predictions for C2H, C3H2 and 
C4H. For both models, the abundance of hydrocarbons 
presents two peaks, the first one located near Ay = 2, 
and a second one near Ay = 5. We shall compare our 
observed data with peak values in the first peak, which 
presents the largest abundances. Given the low spatial 
resolution of the data in comparison with the expected 
spatial variations, the observed abundances are lower lim- 
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02468 10 02468 10 

Ay (mag) Ay (mag) 



Fig. 12. Prediction for the abundances of various species derived from an updated Le Bourlot et al.'s PDR model, using 
two chemistry networks, namely Herbst et al.'s New Standard Model (left) and UMIST95 of Millar and collaborators 
(right). We assumed a gas density of uh = 2 x lO^cm"^ and a radiation field Go = 100. Both models predict too low 
abundances for C-C3H2 and C4H, but also that C2H is present in the UV illuminated gas as observed. 



its to the true abundances. For the Horsehead nebula we 
have good evidence from the maps that emission lines of 
C2H and related molecules peak in the same region as 
the mid-IR emission and that the present data could not 
resolve this sharp peak. 

In both models, the peak abundances relative to H2 
are [C2H] ^ 2 x 10-^ [C3H2] ~ 2 x 10''^° and [C4H] 
~ 2 X 10^1". The peak CO abundance is ~ 2.6 x 10"^, 
since all gas phase carbon is locked into CO for Ay larger 
than 2 magnitudes. Compared to the observed data, mod- 
els predict similar abundances relative to H2 for C2H, and 
CO ; also C4H and C3H2 behave in the same way as we 
observed. However, their predicted abundances relative to 
H2 are lower by an order of magnitude than abundances 
derived from the observations. The discrepancy is clear in 
both models, all the more if we consider that the peak 
abundance is obtained over a fairly narrow range of Ay . 
Since our data are taken with a low spatial resolution, we 
suffer from beam dilution as explained above. How can we 
explain this discrepancy ? We suggest three possibili- 
ties which are not independent: 

— The photo-dissociation rates for hydrocarbons 
(molecules, radicals and ions) included in the chemical 



networks may be overestimated by one order of mag- 
nitude. For most species, the values used in models 
are indeed best guesses with large uncertainties (see 
van Dishoeck 1988 and Fosse 2003). 

— Although including thousands of reactions, chemical 
networks used for models may not list the relevant 
chemical reactions for forming carbon chains and rad- 
icals yet. A related possibility is that some reaction 
rates are incorrect. For example the role of atomic oxy- 
gen in destroying carbon chains and cycles in the in- 
terstellar medium is still not settled (E. Herbst, priv. 
com.). There is also the possibility that neutral 
neutral reactions involving neutral carbon, CH 
, C2H and small hydrocarbons may be at work 
and contribute to the formation of the larger 
hydrocarbons. Reactions of atomic carbon with 
acetylene, ethylene, propyne and allene have in- 
deed been studied by Chastaing et al. (2001) 
and shown to be rapid. 

— Another mechanism for producing carbon chains is 
needed, in addition to gas phase chemistry. An at- 
tractive possibility is the erosion and/or destruction of 
PAHs and small carbon grains. Indeed the main pro- 
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cess constraining the smallest possible size of PAHs 
is the destruction of their carbon skeleton. Various 
mechanisms have been proposed (see Allain et al. 1996, 
Scott et al. 1997, Le Page et al. 2001, 2003, Verstraete 
et al. 2001) which all lead to the production of small 
carbon fragments (C2H2, C2 or C3,...). Once released 
in the gas phase, these fragments may further react 
and take part in the synthesis of the carbon chains 
and radicals we observed. Given the large fraction of 
total carbon in solid form (PAHs can represent up to 
10% of total C) compared with the rather small abun- 
dance of carbon chains (see Tab. the production 
of carbon clusters and molecules from the PAHs is a 
plausible mechanism. Such a mechanism v^ras also 
mentioned by Puente et al. (2003) to explain 
the high hydrocarbon abundances observed in 
NGC7023 and the Orion Bar. 

7. Conclusions 

We have conducted an extensive search for small hydrocar- 
bon chains and rings in three Photon-Dominated Regions. 
The map obtained in the millimetre transitions of C2H, c- 
C3H2 and C4H show that all three hydrocarbons are found 
until the border of the PDRs, despite the strong radiation 
fields (Go ^ 100 — 1000). They moreover correlate well 
with each other, in agreement with trends reported in dif- 
fuse clouds (Lucas & Liszt 2000). In the Horsehead nebula, 
the hydrocarbon large-scale distribution seems closely re- 
lated to the emission of aromatic bands observed in the 
mid-IR, a result in contrast with the CO brightness dis- 
tribution. A firm link between the small hydrocarbons 
and the aromatic band carriers however still requires high- 
resolution maps of this area (see Pety et al. 2003). Heavier 
species have been probed at positions of interest. They are 
mostly detected in the Horsehead nebula, with a notice- 
able detection of CeH around 81.8 GHz. To our knowl- 
edge, this is the first detection of this heavy carbon chain 
reported in a PDR. 

We have derived the hydrocarbon column densities 
and abundances relative to H2 and C2H. On the aver- 
age, the column densities of C^^O and of the most abun- 
dant hydrocarbons are similar in all three PDRs, and com- 
parable to results reported in other sources (e.g. Fuente 
et al. 1993, 2003). They exhibit representative ratios of 
order ~ 25, 1 and 15 for [Ci^Ol/pzH], [c-C3H2]/[C4H] 
and [C2H]/[c-C3H2], respectively. We have compared these 
values to column densities reported in the TMC-1 and 
L134N dark clouds, known as efficient carbon factories. 
Although H2 column densities may be underestimated in 
the dark clouds, abundances relative to H2 are similar 
up to chains and rings of 2-3 carbon atoms. The same 
trend is evidenced in carbon budgets derived for PDRs, 
diffuse clouds and dark clouds conditions. Ratios relative 
to C2H however show lower values for the PDRs in com- 
parisons with TMC-1, a very likely consequence of the 
different chemistries at work in respectively illuminated 
and well-shielded clouds. 



PDR model computations (Le Bourlot et al. 1993) us- 
ing two independent gas-phase chemistry networks suc- 
ceed in reproducing some of the observed abundances 
(Ci«0, C2H) and column density ratios ([Ci^Ol/pzH], 
[c-C3H2]/[C4H]). They however fail in producing the ob- 
served amount of C-C3H2 and C4H by an order of magni- 
tude. This discrepancy may be to due to incorrect dissoci- 
ation rates used in the networks, or to missing gas-phase 
reactions in the considered network. An additional forma- 
tion process of hydrocarbon could also be at work. We 
propose to consider further the role of the AIB carriers in 
the chemistry as their erosion by the UV radiation may 
return small carbon molecules in the gas phase (see e.g. 
Allain et al. 1996, Scott et al. 1997, Le Page et al. 2001, 
2003). New laboratory experiments (Joblin 2002) confirm 
that the photo-dissociation of PAHs can generate various 
kind of carbon clusters whose further dissociation leads 
to the release of small fragments, essentially C3 and C2. 
Further understanding and implementation of this mech- 
anism in ISM models is needed in the future. More obser- 
vational data, especially on heavier molecules, would also 
complete the scarce information on the carbon chemistry 
of PDRs. 

Appendix A: Column density computations 

The various column densities considered in the study 
(Sect. 15. 1|) have been computed under distinct approxi- 
mations, depending of the species and/or the source. We 
describe here the assumptions and details of their compu- 
tation. 

A.l. C^O 

We applied the method described in e.g. Abergel et al. 
(2003), consisting in using the optically thick emission 
lines of ^^CO to infer the kinetic temperatures in all three 
sources. The values obtained at the CO emission peaks 
are indicated in Sect. 15.11 Based on these temperatures, 
we used a LVG approximation and fitted, when applicable, 
the ( J=l-0) and (J=2-l) transitions of C^^O with pairs 
of [n-ii2,N{G^^0)]. When only one transition was available 
(e.g. IC63), an additional assumption was made on the 
volume density upja ■ 

A.2. C2H 

We took advantage from the simultaneous detection of 
4 components of the hyperfine structure of the N=l-0 
transition of C2H and applied the so-called HyperFine 
Structure method (HFS, see e.g.) featured in the CLASS 
reduction software. At each position, the fit returns the to- 
tal transition opacity Ttot, as well as the product [Ji/(Tcx) — 
Jiy{Thg)] XTtot (with Tbg the cosmological background tem- 
perature). The method gives satisfactory results pro- 
vided the signal-to-noise ratio is sufficient, and the opaci- 
ties are not too low (rtot > 0.1). 
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For optically thin or marginally optically thick lines, 
the column density is directly proportional to the in- 
tegrated line intensity. For the excitation temperature 
ranges inferred in the three objects, A^(C2H) is not sensi- 
tive to Tex, and we finally adopted the following expres- 
sions (in cm~^): 



(A.l) 



N{C2il) = (7 ± 1) • 10^=^ X W(C2H)/G(t) 

where G(t) = (1 - e-^')/^', t' = 0.679 x r°-9", is a cor- 
rection factor accounting for non-negligible opacity cor- 
rections (Wyrowski et al. 1999), and r is the opacity of 
the main line at 87.316 GHz. 

A.3. C-C3AV2 

Based on LVG simulations, we found that the C-C3H2 in- 
tensities can be very sensitive to density variations in the 
conditions encountered in our targets (1 — 5 x lO^cm"^, 
Tkin — 20 — 50 K). As an illustration, considering the 
Horsehead case, the density uncertainties derived from the 
C^^O analysis would lead to 7V(c-C3H) variations larger 
of 50%. To solve this problem, we used the good corre- 
lation observed between the brightness temperatures of 
C-C3H2 and C2H (e.g. Fig.O. Since C2H is much less af- 
fected by density variations, and considering the low op- 
tical depth of the molecules, this trend suggests homoge- 
neous physical conditions in the clouds, although the exci- 
tation of each species strongly differs. We thus assumed a 
constant volume density and a representative kinetic tem- 
peratures to apply to the LVG simulations. Since our LVG 
code treats ortho and para forms separatly, we applied to 
the ortho c-C3H2(2i_2-lo,i) transition a standard o/p ra- 
tio of 3. 

A.4. C4H 

The low line intensities suggest very likely optically thin 
lines. Applying the same reasoning as for C-C3H2, we per- 
formed LVG simulations on the C4H(9-8) transition lines. 
Since the model does not account for the hyperfine struc- 
ture of this transition, we built synthetic lines from the 
emission of all observed components. This combination 
uses the expected (at LTE) relative line intensity ratios 
ajy, and writes: 



synth 



N 



where N is the n**^ hyperfine component. In practice, the 
relative intensities of the (N=9-8, J=19/2-17/2, hereafter 
labelled 1) and (N=9-8, J=17/2-15/2, labelled 2) transi- 
tions (consisting both of two blended lines) are 0.53 and 
0.47 respectively, so that W^synth ~ VFi -I- W^2. The outputs 
of the LVG simulations are given in Tab. |21 

A. 5. Other species 

The emissions exhibited by the rarer molecules intro- 
duced in Sect.EK^-CaH, C-C3H, 1-C3.il2, 1-C4il2, CH3C2H, 



CeH and HC3N) have been analysed under LTE assump- 
tions, using the densities and excitation temperatures in- 
ferred from the analyses described previously. In partic- 
ular, for C-C3II and I-C3H we used the Tex derived from 
C2H. For CgH, there are no cross sections available. We 
estimated them from those of IIC3N (Green & Chapman 
1978, - see Cernicharo et al. 1999). Finally, for HC3N, 
we assumed the Tkin values inferred in Sect. lA.ll In case 
of non-detection, upper limits of the column densities are 
given (Tab.EJ. 

A. 6. H2 column densities 

Meaningful comparisons of the observed molecular column 
densities require an accurate estimate of the total II2 col- 
umn densities at each of the positions reported in Tab. El 
We discuss here the methods used to derive as accurate 
as possible values for each of our targets. 

A.6.1. L1688-W 

In this source, we cannot use the extinction map of 
Wilking & Lada (1983) as our positions lie outside their 
coverage. We nevertheless reproduced their approach, con- 
sisting in converting C^^O column density into visual ex- 
tinction, and applied the iV(H2)/ylv calibration by Bohlin 
et al. (1978), which assumes a gas-to-dust ratio of 100. The 
choice of the the Ay/N{C^^O) ratio is not straightforward 
as most reported values have been obtained in cold dark 
clouds where depletion onto grains affects the measure- 
ments and consequently increase this ratio (e.g. Frerking 
et al. 1982, Cernicharo et al. 1987). One can however re- 
compute this ratio for warmer conditions assuming that 
all CO is in the gas phase. Using high-resolution infrared 
spectroscopy of H2 and CO in NGC2024 IRS2, Kulesa et 
al. (2002) recently reported a new measurement of the 
CO abundance relative to H2, [CO] = 3 x 10"^. With 
[Oi^]/[Oi8]=560 (Wilson & Rood 1994) and using stan- 
dard conversion numbers, one can infer: 



Ay/N{C^^O) = i? x 6.6 X 10" 



mag cm 



(A.3) 



where R = Ay/E{B — V) has reported values ranging 
between ~ 3 — 5 (Spitzer 1978, Ducati et al. 2003 and ref. 
therein). This yields: 



(A.2) Ay/N{C^^O) = 2 - 3.3 X 10"^^ magcm^ 



(A.4) 



in very good agreement with the ratio proposed by 
Wilking & Lada (1983). A zero point offset in Ay is how- 
ever to be considered to take into account the fact that 
C^*0 is photo-dissociated for low extinctions. We used 
an offset of 1.5 mag, consistent with the calibrations by 
Frerking et al. (1982) and Cernicharo et al. (1987). We fi- 
nally adopted the following conversion rule for L1688-W: 



iV(H2) (2.5 ± 0.6) • 10" X A^(C'O) + 1.4 • 10^' cm(A.5) 

At the peak C^^O column density, this yields A^(Il2) = 
1.1 X 10^2 cm-2. 
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A. 6. 2. Horsehead nebula 

In the Horsehead nebula, no extinction estimate is avail- 
able at spatial resolution of order 30". With a beam res- 
olution of ~ 2', Kramer et al. (1996) report Ay ~3. This 
is consistent with rough estimates from the IRAS 100 /im 
emission (4' beam) but is very likely not adapted to the 
small-scale structure of this source. Applying the same 
approach as in L1688-W implies extinctions of order 8- 
13 mag. In order to cross-check these estimates, we have 
our observations of the continuum at 1.2 mm (Fig. and 
exploited them in the fashion described by Motte et al. 
(1998), with: 



Si. 



2mm 



,^mHKl.2-Bi.2(rdust) 



(A.6) 



where Smm is the 1.2 mm flux density in mJy per 
11" beam, ribeam the beam solid angle, /i the mean molec- 
ular weight, mn the mass of atomic hydrogen, Bi.2{T) 
the Planck function at wavelength 1.2 mm and tempera- 
ture T, and K1.2 the dust mass opacity at 1.2 mm. Such 
estimates are affected by the approximative knowledge of 
the dust properties at these wavelengths (e.g. Ossenkopf 
& Henning 1994, Stepnik et al. 2003). With the addi- 
tional assumption that Tdust is well represented by Tkin 
(~ 40 here), and adopting K1.2 = 3 — 6 x 10~'^gcm~^, 
we derived from the observed peak flux of 56mJy/beam 
a peak extinction of 12-25 mag. After smoothing to the 
3 mm beam, the corresponding extinction range reduces 
to 9-15 mag, which translate into a column density of 
7V(H2) = 1.1 ± 0.3 X 10^2 cm-2. This is significantly 
above any previous estimates (cf Sect. I2.1|l . but in very 
good agreement with the independent calculation using 
C^^O emission. 

A.6.3. IC63 

In IC63, we simply used the value adopted by Jansen et 
al. (1994) N{B.2) = 5 ± 2 x lO^i cm'^, although it may be 
overestimated at positions away from the CO peak. 
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